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A comprehensive study of the total energy of manganese-rich Heusler compounds using density 
functional theory is presented. Starting from a large set of cubic parent systems, the response 
to tetragonal distortions is studied in detail. We single out the systems that remain cubic from 
those that most likely become tetragonal. The driving force of the tetragonal distortion and its 
effect on the magnetic properties, especially where they deviate from the Slater Pauling rule, as 
well as the trends in the Curie temperatures, are highlighted. By means of partial densities of 
states, the electronic structural changes reveal the microscopic origin of the observed trends. We 
focus our attention on the magnetocrystalline anisotropy and find astonishingly high values for 
tetragonal Heusler compounds containing heavy transition metals accompanied by low magnetic 
moments, which indicates that these materials are promising candidates for spin-transfer torque 
magnetization-switching applications. 


I. INTRODUCTION 

The spintronics community demands materials with 
uniaxial anisotropy for spin transfer torque-random ac¬ 
cess memory (STT-RAM) applications, as well as for 
fundamental skyrmion-related research or in the field 
of magnetic shape-memory alloys [1]. Such materials 
are needed to improve the functionality and applica¬ 
bility of modern devices or device concepts within the 
scope of mass production or proof of practicality. Ad¬ 
dressing this request withing the class of compounds 
with Heusler and Heusler-like structures, the task is ap¬ 
proached by means of relatives of familiar systems [2, 3]. 
These relatives are the family of Mn 2 -based Heusler com¬ 
pounds, the famous pioneering material and ancestor 
of which is Mn 3 Ga [4, 5]. The uniaxial magnetocrys¬ 
talline anisotropy of Mn 3 Ga has been calculated [6] and 
measured [7] several times on different occasions. It is 
thought that anisotropic materials such as these could 
constitute the foundation for magnetic racetrack mem¬ 
ory as proposed by Parkin et a 1. [8]. The indispensable 
perpendicular magnetic anisotropy (PMA) in (ultra-)thin 
structures is best controlled by intrinsic properties rather 
than by shape- or strain-induced anisotropy. As a re¬ 
sult, the perpendicular orientation of magnetization is a 
desired property of the material. Recently, Mn 2 -based 
Heusler systems, yin 2 YZ, were reconsidered as promis¬ 
ing materials. Thus, considerable research has been done 
on related systems. M^NiGa in particular is a well- 
studied material as it is directly related to Ni 2 MnGa, 
which has been the most studied and best understood fer¬ 
romagnetic shape-memory alloy since its discovery [9,10]. 
Mn 2 NiGa, however, is a ferrimagnetic shape-memory al¬ 
loy that is theoretically linked to N^MnGa through a 
substitution series, with a transition from ferro- to ferri¬ 
magnetic ordering due to the increasing manganese con¬ 
tent. In addition to Mi^NiGa, other Mn 2 -based Heusler 
alloys have been synthesized and characterized or have 
been theoretically treated, such as Mn 2 CoGa [11-15], 
Mn 2 FeGa [16], and Mn 3 Ga [4, 7]. In addition to ma¬ 


terials in which Y is an atom from period IV (or the 
3d series), equivalent systems with heavier constituents 
as Y species have been investigated. Among these were 
Mn 2 RuGa [17] (which has been found to have more or 
less random occupation of sites), M^RhGa [18] (cubic, 
disordered), Mn 2 PtGa [19, 20], and M^Ptln [21], in the 
context of large exchange bias effects. 

Although a detailed study of a single material is an im¬ 
portant task that results in valuable knowledge, the in¬ 
clusion of cluster knowledge into a general concept cre¬ 
ates comprehensive insight. The same intention that 
guided our previous work [15] motivated us to undertake 
a similar approach in the current study, i.e., comprehen¬ 
sively studying a selected set of systems. We intend to 
understand the general trends governing the formation 
and magnetism of tetragonal materials for the aforemen¬ 
tioned applications. 

In this paper, we show how the magnetism and 
the atomic structure change withinthe M^Y^^Ga, 
Mn 2 Y (4d) Ga, and Mn 2 Y (5d) Ga series. The trigger quan¬ 
tities causing the tetragonal distortion, as well as the con¬ 
sequences of this distortion, i.e., the magnetocrystalline 
anisotropy (MCA) [3] , will be highlighted and placed into 
an appropriate context. 


II. CRYSTAL STRUCTURE 

Heusler alloys are nowadays informally divided into 
two structure types, the so-called regular and in¬ 
verse types, referring to the original Heusler compound 
Cu 2 MnAl as the reference system [22]. The materials 
first associated with the name of Fritz Heusler were cubic 
phases analogous to Cu 2 MnAl, with the stoichiometry 
X 2 YZ. From then on, similar materials were thus labeled 
“Heusler compound”, extending the original definition of 
Heusler compounds to the family of Heusler materials 
that incorporates a variety of similar structures. These 
structures are derived from the original compound, with 
occupation of the Wyckoff positions, 8c, 4b, and 4a, in 
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FIG. 1. (Color online) The conversion from a (a) cubic Heusler structure to a (b) tetragonal derivative phase is displayed in 
terms of a non-displasive transformation for a system with the general composition of XX'YZ (X, X', Y - transition metals 
and Z - main-group element, marked as red, orange, blue, and green, respectively) within the fee lattice, (c) Relationship of 
the lattice parameters as atet = a cu b/\/2. 


the spacegroup (SG) 225 through the introduction of va¬ 
cancies or slight structural changes. These modifications 
alter the structure by breaking the inversion symmetry 
when going from a regular Heusler phase to an inverse or 
/iaZ/-Heusler structure. 


Heusler 

SG 

4 d 

4c 

4 b 

4a 

regular 

L2i 

225 

Mn 

Mn 

Y 

Ga 

inverse 

x a 

216 

Mn 

Y 

Mn 

Ga 

half 

Clb 

216 

Mn 

Mn 

□ 

Ga 


TABLE I. The structural relationship of the regular and 
inverse structure types for Mn 2 VGa is shown, where 
X=X'=Mn and Z= Ga, as compared to Fig. 1. Only one 
possible configuration of a ftaZ/-Heusler type is listed exem¬ 
plary. 

Heusler materials are generally understood to be in- 
termetallic compounds, distinguishing them from gen¬ 
eral intermetallics which are forming a broad range of 
solid solutions, with no preferred but statistical occupa¬ 
tion of crystallographic sites. They are also set apart 
from other ionic or covalent compounds because Heusler 
systems allow the formation of substitutional series of 
single sites. However, they maintain the character of an 
ordered compound, and thus they are on the borderline 
between alloys and compounds. Some distinct systems 
exhibit a tendency to form alloys, nevertheless. The gen¬ 
eral composition is given by XX'YZ , where the classical 
definition has been widened to incorporate quaternary 
materials within the family of Heusler compounds. In the 
representation of SG 216 (E43m), the structure contains 
four highly symmetric Wyckoff positions: 4 d (§,§,§), 

4c ib b l)’ Ab (|> b D> and 4o (°>°>°)- Depending on 
the occupation of the crystallographic positions, two or¬ 
dering possibilities for ternary alloys ( X = X') are ob¬ 
tained (Table I). In this study, gallium was chosen as 
the Z element, whereas one manganese atom, X , occu¬ 
pies position 4 d. The second manganese atom, X\ and 
the other transition metal, Y , are located at 4c or 4&, 


respectively, depending on the formation of the regular- 
or inverse-type Heusler material, as shown in Table I. 
Other Heusler-related structures [23] are the tetragonal 
derivatives of the cubic parent phases, which have been 
widely treated in the context of magnetic shape-memory 
alloys. The relationship and the unit cell transformation 
between cubic and tetragonal phases is depicted in Fig. 1. 
It is seen that a conventional cubic unit cell can be de¬ 
scribed in terms of a tetragonal lattice exhibiting a c/a 
ratio of \/2. The cell parameters are interrelated accord¬ 
ing to Ctet = Ccub , atet = a cub /V2- In this study, a set of 
Mn 2 -based materials, including transition metals of peri¬ 
ods IV, V, and VI, are considered, namely, Mn 2 V^ 3d ^Ga, 
Mn 2 V (4 <i)Ga, and Mn 2 F (5[i) Ga. 


III. COMPUTATIONAL DETAILS 

The numerical work was done within density func¬ 
tional theory as implemented in the all-electron FP- 
LAPW code, WIEN2K [24], employing the generalized 
gradient approximation (GGA) in the parametrization 
of Perdew, Burke, and Enzerhof as exchange-correlation 
functional [25]. The angular momentum truncation was 
set to Z max = 9 and the number of plane-waves deter¬ 
mined by RK max = 9 to ensure well-converged calcu¬ 
lations. The energy convergence criterion for the self- 
consistent field calculations was set to 10 -5 Ry, whereas 
the charge convergence was set to 10 -3 . All calculations 
where done on a 20 x 20 x 20 fc-mesh. For a set of given 
c/a ratios, the volumes were optimized and fitted to the 
Birch-Murnaghan equation of state [26, 27]. From these, 
the optimal ratio and volume were obtained and the lat¬ 
tice parameters evaluated. On the basis of the crystallo¬ 
graphic details, complex magnetic properties such as the 
exchange parameters and the corresponding Curie tem¬ 
peratures (Tq ) were computed by means of the Korringa- 
Kohn-Rostoker (KKR) Green’s function method as im¬ 
plemented by the Munich SPR-KKR package [28]. The 
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FIG. 2. (Color online) Energetic response to volume-conserving elongations and compressions of the crystal structure along 
the c-axis for the M^F^G a-series. Panel (a) includes systems with a valence electron count Ny (Y) < 7 (early transition 
metals, ETM) and panel (b) shows the remaining combinations of M^F^^Ga with Ny (F) > 7, with so-called late transition 
metals (LTM). 


angular momentum expansion was truncated for Z max = 3 
which corresponds to the /-wave symmetry. The energy 
integration was done on a complex energy mesh with 48 
points along the integration path, using the Lloyd’s for¬ 
mula [29] for an improved estimate of the Fermi energy. 
The computation of the exchange parameters is based 
on the classical Heisenberg model, which was evaluated 
by means of the real-space approach [30] . This provides 
site- and distance-dependent exchange between sites via 
infinitesimal rotation of the magnetic moments at a par¬ 
ticular site in real space. To account for distance depen¬ 
dence, an appropriate truncation of the cluster radius, r, 
around each atomic site had been chosen. This radius 
was set to 3.5 a lattice spacings to capture even small 
interactions, as the largest contributions to the effective 
exchange constants are found for radii smaller than 1.5 a 
lattice spacings [15, 31]. 


IV. RESULTS 

A. Lattice Relaxation 

The total energy E (c/a) as a function of the c/a ratio 
was calculated and the results are shown in Figs. 2-5. 
The energy zero are defined with respect to the cubic 
parent compound, and consequently the energy differ¬ 
ence for all phases can be compared easily. The case of 
Mn 3 Ga is used as a benchmark and is repeatedly plotted 
in Figs. 2 and 3. Table III contains the numerically opti¬ 
mized lattice parameters. The study reveals that a large 
number of the herein-treated materials are most stable 
in their respective tetragonal structures, with c/a > y/2. 


Tetragonally compressed structures are described by 
c/a < v/2, whereas tetragonally elongated lattices are 
characterized by c/a > y/2, as compared to the cubic 
parent or austenite phase. Elongation occurs with an 
increase in the length of the c-axis, whereas the ab- 
plane is compressed, leaving the volume approximately 
unchanged. In this study, the volume of the unit cell, 
14eii, was optimized in addition to the c/a- ratio, and it 
was found that no significant change occurred for most 
cases (Table II). 

Figs. 2a and 2b show the E(c/a) curves of the 


Material 

C t/ c 

Material 

Ct/c 

Mn2ScGa 

5.08 

Mn 2 OsGa 

2.32 

Mn 2 CrGa 

6.11 

Mn 2 IrGa 

2.44 

Mn 3 Ga 

9.08 

Mn2PtGa 

-3.13 

Mn 2 FeGa 

1.93 

Mn2AuGa 

-1.49 

Mn 2 CoGa 

1.82 

Mn 2 OsSn 


Mii2NiGa 

-1.09 

Mn 2 IrSn 

-1.94 

M 112 CuGa 

2.37 

Mn 2 PtSn 

-0.35 

Mn2MoGa 

3.46 

Mn20sln 

-0.36 

Mn2RuGa 

1.63 

Mn2lrln 

-1.01 

Mn2RhGa 

1.48 

Mn2PtIn 

-0.86 

Mii2PdGa 

-2.78 

Mn 3 Ge 

3.98 



Mn 2 FeGe 

3.80 


TABLE II. Relative volume change between the austenite and 
martensite phase in percent, C t / c = (Vtet — Kub) • lOO/Kub 

Mn 2 Y (3d) Ga series for early transition metals (ETM, 
Y = Sc, Ti, V, Cr) and late transition metals (LTM, 
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FIG. 3. (Color online) Energetic response to volume-conserving elongations and compressions of the crystal structure along 
the c-axis for the (a) M^Y^^Ga and (b) M^Y^^Ga series. Systems involving late transition metals and systems exhibiting 
at least a local tetragonal energy minimum are shown. 


Y = Mn, Fe, Co, Ni, Cu), respectively. A preference 
for tetragonal structures is seen for materials including 
LTM for M^F^-’Ga, with Y being Mn, Fe, or Ni. 
In the group of ETM systems, M^ScGa was found to 
be tetragonal, which is an exception. The situation for 
M^F^-’Ga and Mn 2 F^ 5d ^Ga is seen to be similar. 

In contrast to M^CoGa, for which a cubic structure 
(c/a = -\/2) is preferred, the systems with the same va¬ 
lence electron count (Ny = 26) involving heavier species 
(ruthenium and osmium) exhibit a global energy mini¬ 
mum for the tetragonal structure. 

Discontinuities in the E (c/a) curves are observed for the 
cases of MnaGa and M^FeGa, resembling a first-order 
transition, whereas E(c/a ) for M^NiGa is continuous 
(Fig. 2b). Strongly composition-dependent modulated 
martensitic phases [32, 33] and premartensitic phases 
[34, 35] have been experimentally produced; in combi¬ 
nation with first-principles calculations [36] , the onset of 
the martensitic transition has been thought to be initi¬ 
ated by a displacement of atomic planes orthogonal to the 
crystallographic c-axis of the tetragonal cell. Experimen¬ 
tally [35] and theoretically [36], it has been shown that 
phonon softening along [£, /, 0] exists in shape-memory 
materials, and thus the transition has been related to 
the occurrence of the tetragonal distortion. 

The approach undertaken in the current study resulted 
in elongated structures in those cases, which were unsta¬ 
ble towards a tetragonal distortion. Additionally, global 
energy minima for distorted phases were found mostly for 
the inverse structures, as depicted in Fig. 4. In contrast, 
compressed variants possessing a global energy minimum 
were not observed in the Mn-Ga Heusler family. The cal¬ 
culated lattice data and the figures shown (Figs. 2-3) lead 
to the conclusion that stable tetragonal structures may 


only be formed in the series of M^YGa that includes 
LTMs, making them derivatives of inverse Heusler sys¬ 
tems. The information obtained from the lattice opti¬ 
mization in terms of the relative positions of the energy 
minima is condensed in Fig. 4, which gives the preferred 
crystal structure visualized in the manner of the periodic 
table. 

Each compound in the family of M^YGa materials 
shown in this figure is depicted by one cell that symbol¬ 
izes a transition metal Y of the 3d, 4d, or 5 d series. The 
corresponding cell is built up by three subcells, which 
represent the two variants of chemical coordination (the 
first or second subccll) and the existence of a global en¬ 
ergy minimum for c/a ^ y/2 (third subcell). The color 
code symbolizes the energy level of a configuration on 
the energy landscape relative to one another: dark green 
- global minimum, light green - local minimum, red - no 
minimum. In cases in which the investigated materials 
did not exhibit a cubic minimum, the first two subcells 
are understood to be the type of coordination around 
the 4 d crystallographic site: the symmetry of the coordi¬ 
nating shell is either centrosymmetric (in relation to SG 
225) or non-centrosymmetric (in relation to SG 216). As 
shown in a preceding publication on the cubic variants 
[15], systems involving ETMs adopt the L2i-type struc¬ 
ture, whereas compounds containing LTMs are found to 
have the inverse Heusler structure (X a -type). 

Inspection of Fig. 4 in combination with Figs. 2 and 3 
reveals interesting details, such as the fact that tetrag¬ 
onal derivative phases of cubic Heusler alloys, implying 
that a global energy minimum is present, are observed 
only for a valence electron count of Ny > 24. It is also 
clearly seen that the onset of the formation of tetrago- 
nally elongated structures evolves over the periods from 
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FIG. 4. (Color online) Schematic overview of the preferred site occupancy and crystal structure of Mi^TGa Heusler compounds. 
Stable, metastable, and instable lattices are marked by dark-green, light-green, and red subcells, respectively. 


left to right and from lower to higher Ny- 
This leads to the question of the mechanism behind this 
distortion, which we are going to approach in Sec. IV A1. 
The lattice parameters for all cubic compounds are found 
within a range of Aa c /^ = 0-35 A for the Ga series. The 
values increase from the borders of the series towards 
the middle of the range. The same behavior is found 
for tetragonal compounds, where the range spans from 
atet = 3.68 to 4.11 (Aatet = 0.43 A) and c te t = 6.91 to 
7.48 (Ac te t = 0.54 A), whereas Mr^ScGa is an exception. 
The c/a coordinate exhibits inverted behavior, decreas¬ 
ing from the middle to the left and right borders of the 
series. The overall similarity of the lattice data opens the 
possibility of intermixability with each other and thus 
tunability of the whole class of materials. Therefore, the 
magnetization may be adjusted over a large range, al¬ 
lowing for the formation of tetragonal compensated fer- 
rimagnets via adequate substitution. 

Having discussed the structural trends as a function 
of the valence electron change in the d-electron sys¬ 
tem, the effects of the variation of the main group el¬ 
ement, Z , from Ga to Ge, In, and Sn in a small sam¬ 
ple of compounds [Mn 2 (Mn,Fe)Ge, Mn 2 (Os,Ir,Pt)In, and 
Mn 2 (Os,Ir,Pt)Sn] are briefly discussed. The correspond¬ 
ing data, including those for MnaGa, are graphed in 
Fig. 5. 

Fig. 5 reveals that the placement of In and Sn at the Z po¬ 
sition leads to the emergence of a deep energy minimum 
for Mn 2 lrSn with c/a < \[2 and even deeper minima 
for elongated phases with large c/a-ratios of Mn 2 OsIn, 
Mn 2 IrIn and Mn 2 OsSn whose lattice parameters resem¬ 
ble those of layered structures. Heusler alloys are of¬ 
ten interpreted in terms of a rigid, band model-like ap¬ 
proach, and an electron-filling scheme is employed for 
the prediction of magnetic moments in Co 2 -based al¬ 
loys. Interpretation of the tetragonal instabilities us¬ 
ing such an approach leads to the assumption of MnaGe 


(Ny = 25) behaving analogously to Mn 2 FeGa (Ny = 25) 
and of Mn 2 FeGe (Ny = 26) being similar to Mn 2 CoGa 
(Ny = 26). A comparison of the corresponding c/a 
curves demonstrates that this is approximately true and 
thus that isoelectronicity is an appropriate concept in the 
chemistry and physics of Heusler materials. 

Although our calculations agree well with those of previ¬ 
ous works [37], comparison with the experimental data of 
Mn-Ni-Ga systems exhibits an interesting discrepancy. 
The deviation between experiment and theory was traced 
back to the deviating structural model in terms of the oc¬ 
cupation of the involved sites. Neutron diffraction studies 
on Mn 2 NiGa highlighted the fact that the order is differ¬ 
ent from the expected X a -type in the austenite phase. 
Thus, the chemical formula reads as (Mn,Ni) 2 MnGa and 
is called L2if,-type because the point symmetry includes 
inversion symmetry through random occupation of 4 d 
and 4c sites with Ni and Mn [38]. Similarly in Mr^FeGa, 
a deviation from theory was found in an experimental 
study because the site occupation was expected to differ 
from the perfect MnFeMnGa ordering [16]. Similar issues 
are present for the M^H^^Ga series. From available 
data, including neutron diffraction studies, the site occu¬ 
pancies have been clarified. Several authors found that 
members of the M^RuZ series (Z = Ga, Ge, Si, Sn) 
and Mn 2 RhZ series (Z = Ga) [17, 18, 39] not to be 
tetragonal under the respective reaction conditions. In 
contrast, they have been realized as cubic Heusler alloys 
exhibiting a strong degree of anti-site disorder, which has 
been characterized as an alloying tendency [18]. 
Orthorhombic deformations of the unit cell have been 
observed in some systems such as Mn 2 NiGa [37]. Thus, 
the restriction to tetragonal distortions and ordered com¬ 
pounds in conducted studies leads to a simplified descrip¬ 
tion of these materials. Nevertheless, a general under¬ 
standing can be obtained in approaching the Mn 2 -based 
Heusler systems through this ansatz. In future studies 
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that aim to predict ground-state structures and magnetic 
configurations, the parameter space for the atomic sites 
and relaxations paths for the electronic and spin degrees 
of freedom has to be enlarged and restrictions that are 
widely used have to be dropped. As a consequence of 
the applied restrictions, disorder effects were not incorpo¬ 
rated into this study, thus leaving open any explanations 
of the deviations from experimental results. 


c/a (cubic reference cell) 

0.8 1.0 1.2 1.4 1.6 1.8 



FIG. 5. (Color online) The energy landscapes, E(c/a), 
of MnaGe, Mii 2 FeGe, Mii 2 (Ir,Pt)Sn, and Mn 2 (Ir,Pt)In are 
shown. 


1. Analysis of the Densities of States 

Various attempts to explain the instability of the cu¬ 
bic phase have been given in the literature using mod¬ 
els such as the band Jahn-Tcller effect [5, 16], anoma¬ 
lous phonon vibrations [49, 50], and Fermi-surface nest¬ 
ing [51]. These different approaches describe the same 
behavior, i.e., the instability of the cubic phase, from 
different perspectives and extracting different types of 
information. Commonly, the densities of states (DOSs) 
of related austenite and martensite are compared and 
contrasted marking the starting and the endpoint of the 
transition. As example for the Mn 2 Y^ 3d ^Z series, the 
DOSs of MnaGa, M^FeGa, and Mi^NiGa are shown 
in Figs. 6-8. For the cubic variants, the partial DOSs 
(PDOSs) are shown in their corresponding projections on 
the sites and in terms of projections on the irreducible 
representations. The peaks in the PDOS of Mn(46) are 
well separated on the energy scale. On the one hand, 
this separation is due to the strong crystal field split¬ 
ting of the Mn(45) d-states, where the occupied e g states 
are located in a range between —4 and —3 eV, whereas 
the t 2 g PDOS is found between —1 eV and the Fermi 
edge, £f ■ On the other hand, the separation of occu¬ 
pied and empty states follows from the strong exchange 
split of the Mn(46) d-states, in contrast to the d-states 


of Y (4c) and Mn(4d), which are found to be more widely 
dispersed even though the majority and minority states 
are separated owing to exchange splitting. The PDOSs 
of Mn(4d) and Y (4c) are strongly dispersed, with con¬ 
siderable overlap of the spectral weight between the t 2g 
and e g states in the majority channel, whereas the minor¬ 
ity channel is gapped, with t 2g characterizing the lower 
boundary and e g comprising the upper boundary of the 
gap. As the DOS is gapped in the minority spin channel 
the study of the tetragonal system is strongly facilitated 
as the states at the Fermi edge in the majority channel 
mostly constitute the origin of the tetragonal distortion. 
The majority PDOSs in the range of —5 eV up to £f 
exhibit a characteristically shaped peak structure. From 
Mi^Ga to Mn 2 CoGa, the majority spin channel (lower 
panels of Figs. 6a-8c) is continuously filled. The Fermi 
energy is consequently shifted to higher band energies 
and thus Ef sweeps over a range of the majority DOSs, 
whereas the minority spin channel remains unchanged. It 
is clearly seen from Figs. 6 and 7 that the tetragonal tran¬ 
sition correlates with the peak structure of the majority 
DOSs. For Ep being centered on a peak of the majority 
DOS, the tetragonal distortion can be triggered. These 
local maxima are mainly composed of states of the Y (4c) 
(Mn, Fe, Co, Ni) and Mil (4 d) atoms. These energy lev¬ 
els are of t 2g symmetry mainly. In simple interpretation, 
the DOS can be understood in a rigid-band-like fashion. 
The limit of this interpretation is reached with Mi^NiGa, 
where the d-PDOS is rearranged and the Slater-Pauling 
rule is no longer valid for the cubic phase [15]. Compar¬ 
ing this to the PDOS of the tetragonally distorted sys¬ 
tems, it is observed that the resulting PDOSs are widely 
dispersed and significantly less structured. Mn 2 NiGa, 
however, behaves differently. Further filling of the ma¬ 
jority states, as intuitively expected, does not occur. In 
contrast, the gap in the minority channel closes as the 
states of e g symmetry are pulled towards ef ■ Thus, the 
tetragonal distortion in Mn 2 NiGa is formed by another 
mechanism, which may explain why M^NiGa is found 
to be a magnetic shape memory alloy, whereas M^Ga 
and Mn 2 FeGa are found in their respective tetragonal 
crystal structures, although the total energy differences 
are comparable to that of M^NiGa. 

We emphasize that these findings differ from other mod¬ 
els in which the instability is thought to depend solely 
on states of Mn(46) that is found in a tetrahedral envi¬ 
ronment. The DOS at Ef in mainly composed of states 
of the 4c and 4 d positions with minor contributions from 
Mn(46). Here the instability removes the strong peaks of 
Mn, Fe, and Co at the 4c position, whereas the states of 
Mn(46) are not rearranged. 


2. Spin-Polarization 

Half-metallicity [52, 53] (complete or nearly complete 
spin-polarization P(sf) ~ 100%) is generally observed 
in cubic C 02 - and Mn 2 -based Heusler compounds. The 
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TABLE III. Calculated lattice parameters of the cubic parent phases and corresponding tetragonal structures, in comparison 
with existing literature data. The lattice parameters are given in Angstroms. The saturation magnetization data, Ms, are 
given in pe/ f.u. The listed literature data refer to theoretical (xt) or experimental (xe) investigation of the mentioned cubic 
( cy ) or tetragonal (t y) phases, respectively. 


N v 

SG 

c/a a c 

a t 

M s 

SG 

c/a 

at 

Ct 

M s 

SG 

c/a 

a t 

Ct 


ordering 

Lit. 

Mn2ScGa 

20 

225 

s/2 6.15 

4.35 

-4.00 

139 

1.94 

3.98 

7.70 

-5.36 








Mn2TiGa 

21 

225 

s/2 5.95 

4.21 

-2.97 

139 





225 

s/2 


5.95 

ct 

MnMnTiGa 

[40] 

Mn 2 VGa 

22 

225 

s/2 5.82 

4.12 

-1.98 

139 





225 

s/2 


5.91 

ce 

MnMnVGa 

[41],[42] 

Mn 2 CrGa 

23 

225 

s/2 5.76 

4.07 

-1.00 

139 

1.82 

3.82 

6.95 

-2.75 

225 

V2 


5.77 

ct 

MnMnCrGa 

[43] 

MnsGa 

24 

225 

s/2 5.82 

4.12 

0.01 

139 

1.82 

3.90 

7.08 

-1.89 

139 

1.77 

3.77 

7.16 

tt 

MnMnMnGa 

[4] 












139 

1.82 

3.90 

7.09 

te 

MnMnMnGa 

[4] 

Mn 2 FeGa 

25 

216 

s/2 5.79 

4.09 

1.03 

119 

1.98 

3.68 

7.29 

-0.78 

119 

1.90 

3.79 

7.19 

te 

Mn(Fe,Mn) 2 Ga 

[16] 

Mn2CoGa 

26 

216 

s/2 5.78 

4.09 

2.00 

119 

1.93 

3.71 

7.14 

0.17 

216 

s/2 


5.86 

ce 

MnCoMnGa 

[44] 

Mn2NiGa 

27 

216 

s/2 5.85 

4.14 

1.18 

119 

1.82 

3.79 

6.91 

1.00 

216 

s/2 


5.91 

ce 


[45] 












119 

1.72 

3.91 

6.70 

te 


[45] 












225 

s/2 


5.94 

ce 

(Mn,Ni) 2 MnGa 

[38] 












139 

1.74 

3.89 

6.77 

te 

(Mn,Ni) 2 MnGa 

[38] 












216 

s/2 


5.84 

ct 

MnNiMnGa 

[37] 

Mn2CuGa 

28 

216 

s/2 5.94 

4.20 

0.33 






216 

s/2 


5.94 

ct 

MnCuMnGa 

[37] 

Mn2ZrGa 

21 

225 

s/2 6.14 

4.34 

-3.00 

119 












Mn 2 NbGa 

22 

225 

s/2 6.00 

4.24 

-2.00 

119 












Mn2MoGa 

23 

225 

s/2 5.91 

4.18 

-1.01 

119 

1.81 

3.89 

7.04 

-2.99 








Mn2RuGa 

25 

216 

s/2 5.96 

4.22 

1.03 

119 

1.96 

3.80 

7.45 

-0.24 

216 

s/2 


6.00 

ce 

(Mn 2 ,Rui )sGa 

[17] 

Mn2RhGa 

26 

216 

s/2 5.98 

4.23 

1.64 

119 

1.94 

3.82 

7.43 

0.10 

225 

V2 


6.03 

ce 

(Mn,Rh) 2 MnGa 

[18] 

Mn2PdGa 

27 

216 

s/2 6.12 

4.33 

0.55 

119 

1.84 

3.93 

7.23 

0.93 








Mn2AgGa 

28 

216 

s/2 6.22 

4.40 

0.34 

119 












Mn 2 HfGa 

21 

225 

s/2 6.12 

4.33 

-2.99 

119 












Mn2TaGa 

22 

225 

s/2 6.00 

4.24 

-1.99 

119 












MnaWGa 

23 

225 

s/2 5.92 

4.19 

-0.94 

119 












Mn20sGa 

25 

216 

s/2 5.95 

4.21 

1.02 

119 

1.97 

3.80 

7.48 

-0.28 








Mn2lrGa 

26 

216 

s/2 5.97 

4.22 

2.00 

119 

1.95 

3.83 

7.44 

0.11 








Mn2PtGa 

27 

216 

s/2 6.13 

4.33 

0.44 

119 

1.87 

3.91 

7.31 

0.75 

119 

1.38 

4.37 

6.05 

te 

MnPtMnGa 

[19] 

Mn2AuGa 

28 

216 

s/2 6.26 

4.42 

0.19 

119 

1.73 

4.11 

7.13 

0.14 








Mn20sSn 

25 

216 

s/2 6.21 

4.39 

1.50 

119 

1.95 

3.97 

7.75 

-0.02 








Mn2lrSn 

26 

216 

s/2 6.31 

4.46 

0.41 

119 

1.91 

4.01 

7.67 

0.45 

119 

1.54 

4.29 

6.59 

te 

MnlrMnSn 

[46] 

Mn2PtSn 

27 

216 

s/2 6.39 

4.52 

0.19 

119 

1.81 

4.15 

7.52 

-0.02 

119 

1.35 

4.51 

6.08 

te 

MnPtMnSn 

[1] 

Mn20sln 

25 

216 

s/2 6.26 

4.43 

0.62 

119 

2.02 

3.93 

7.93 

-0.27 








Mn2lrln 

26 

216 

s/2 6.30 

4.45 

0.68 

119 

1.98 

3.97 

7.85 

0.07 








Mn2PtIn 

27 

216 

s/2 6.37 

4.51 

0.31 

119 

1.84 

4.12 

7.57 

0.38 

119 

1.57 

4.32 

6.77 

te 

MnPtMnln 

[21] 

MnsGe 

25 

216 

s/2 5.76 

4.07 

1.01 

119 

1.90 

3.74 

7.10 

-0.98 

225 

1.91 

3.81 

7.26 

te 

MnMnMnGe 

[47] 












225 

1.90 

3.75 

7.12 

tt 

MnMnMnGe 

[6] 

Mn2FeGe 

26 

216 

s/2 5.73 

4.05 

2.01 

119 

2.05 

3.63 

7.42 

-0.06 

216 

s/2 


5.80 

ct 

MnFeMnFe 

[48] 


highly symmetric structure and the peculiar electronic of a pseudo -gap in one spin-channel is observed instead, 

properties due to covalent bonding lead to the appear¬ 
ance of a gap in the minority density of states (DOS). 

The emergence of the tetragonal distortion reduces the 
spin-polarization of the half-metallic cubic parent phases. 

The degeneracy of the t 2g and e g states is lifted due to 
the change in local coordination caused by the distortion, 
which is seen in Figs. 6-8 and Fig. 16, and the emergence 
















E- e F [eV] 


FIG. 6. (Color online) DOS of cubic (SG 225) and tetragonal 
(SG 139) Mn 3 Ga. 


B. Magnetic Ground State 

1. The Slater-Pauling Rule 

Referring to a previous publication [15] on the cu¬ 
bic parent compounds of the investigated materials, the 
results are presented by means of the Slater-Pauling 
curves. As can be seen in Fig. 10, the Slater-Pauling 
rule experiences strong changes, so that calling it the 
Slater-Pauling rule is done for reasons of convenience. 
Fig. 10 visualizes the fact that the magnetization of 
all tetragonal alloys of the Mn 2 yG)Ga family experi¬ 
ences a shift to smaller values. In the case of the 
Mn 2 y( 3d )Ga group, this shift is found to be constant 
throughout the set of compounds, which results in a 
linear dependence of the net moment on the valence 
electron count, thus giving rise to pseudo-Slater-Pauling 
behavior, even though half-metallicity and thus integer 
net moments are not observed (Fig. 9). A decrease of 
the net moment is also found for the Mn 2 F( 4d )Ga and 
Mn 2 F*- 5d - ) Ga compounds. Unlike the lighter compounds, 
the changes are not constant over the series and therefore 



-7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 

E- s F [eV] 


FIG. 7. (Color online) DOS of cubic (SG 216) and tetragonal 
(SG 119) MnaFeGa. 


result in nearly vanishing net moments for Mn 2 RuGa, 
Mn 2 RhGa, Mn 2 PdGa, Mn 2 OsGa, Mn 2 IrGa, Mn 2 PtGa, 
and Mn 2 AuGa, as is seen in Fig. 10. Nevertheless, com¬ 
pensation of spin moments may be achieved for an elec¬ 
tron count close to Ny = 25.7, which can be realized by 
intermixtures of stoichiometric phases of Mn 2 FGa, such 
as TOMn 3 Ga = —1-89 p# and a corresponding proportion 
of Mn 2 NiGa with 77iMn 2 NiGa = 1-00 pb or Mn 2 CoGa. 
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E- s F [eV] 


FIG. 8. (Color online) DOS of cubic (SG 216) and tetragonal 
(SG 119) Mn 2 NiGa. 


Thus, fractions of 0.435 WMn 2 NiGa and 0.565 mMn 2 FeGa 
could ideally lead to complete compensation of the mag¬ 
netization. A similar approach was undertaken by Nayak 
et a 1. that obtained a compensated ferrimagnet by vary¬ 
ing the Mn/Pt-ratio in Mn 3 _ x Pt x Ga leading to com¬ 
plete compensation of magnetization for x « 0.59 theo¬ 
retically [20]. In Mn 2 Ru a; Ga thin films the compensation 
of the spin-moment has been achieved through variation 


Group Number M 


3456789 10 11 



FIG. 9. (Color online) The spin-polarization at the Fermi en¬ 
ergy for both cubic and tetragonal phases. For a system with 
A r v < 27, the spin-polarization is reduced by the tetragonal 
distortion. 


16 20 24 28 16 20 24 28 32 



FIG. 10. (Color online) Slater-Pauling curves of 
(a) Mn 2 y^ 3l ^Ga, (b) M^T^-’Ga compounds, and (c) 
Mr^T^Ga compounds. 


of the ruthenium concentration [54]. 
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(a) Mn 2 K (3d) Ga (b) Mn/ (4d| Ga (c) Mn 2 K (5d) Ga 

Sc Ti V Cr Mn Fe Co Ni Cu Mo Tc Ru Rh Pd Ag Os Ir Pt Au 



FIG. 11. (Color online) Atomic magnetic moments in (a) Mn 2 >^ 3< ^Ga, (b) Mi^V^Ga, and (c) Mn 2 l^ 5t ^Ga compounds. 
Open symbols denote the L2i-type coordination. Filled symbols denote the X a -type coordination. The squares stand for the 
cubic materials, whereas the tetragonal systems are represented by rectangles. 


2. Local Magnetic Moments 

The change of the total magnetic moment is to be un¬ 
derstood in terms of the site moments. Inspection of 
Fig. 11 clearly reveals the dependencies. The change in 
the local moments indicates the specific impact of the dis¬ 
tortion on the single sites. The local moments responding 
strongest to elongation/compression of the crystal axes 
are members of the magnetic sublattices formed by the 
Mn(8c), Mn(4d), and V(4c) sites (the former so-called 
tetrahedral sites as compared to SG 225). Depending 
on the characters of the local moments, which are either 
itinerant or localized in nature [12, 52], we find major dif¬ 
ferences in the influence of the tetragonal distortion on 
these sites. In Fig. 11, the localized moment of Mn(46) is 
found in the upper part of the plot for a positive value of 
the magnetic moment of approximately 3 fiB ■ The effect 
of elongation along the c-axis and compression of the ab- 
plane has a stronger influence on moments of itinerant 
character found in the lower part of the plot, referring 
to the 4 d and 4c sites. These are located in the same 
lattice plane (compare the atoms depicted as red and 
blue spheres in Fig. 1). Manganese on site 4 b [Mn(46)] 
exhibits a large localized moment of 3 /ib, is thus is gener¬ 
ally much less affected. Apart from the magnitude of the 
local moments, Fe(4c) in M^FeGa exhibits a spin-flip 
from parallel to anti-parallel alignment of the Mn(45)- 
Fe(4c) interaction upon the tetragonal distortion. Apart 
from the changes in magnitude of the local atomic mo¬ 
ments, the effective anti-parallel coupling of the nearest 
neighbor manganese atoms does not suffer from the struc¬ 
tural transformation, as will be quantified in terms of the 
exchange interaction constants in Sec. IV C. Since the 
volume change C t / c is in the order of 1-2% for the most 


systems, the nearest-neighbor distance basically does not 
change, whereas the direction of the nearest-neighbor in¬ 
teraction does (Fig. 16). In the case of Mn 3 Ga, C t / c 
is approximately 9%, and thus may suppress the shape 
memory effect in this material. 


3. Magnetocrystalline Anisotropy 

Inherent in noncubic crystals is a directional preference 
of magnetization that is absent in cubic materials, which 
is related to the tetragonal modification of the crystal 
axes. The magnetocrystalline anisotropy (MCA) energy 
is here defined as the energy difference between states 
with magnetization pointing along the z-axis and the 
x- or y- axis, that is, Emca = -E'(ioo) — -E(ooi)i whereas 
other crystallographic directions are not considered. The 
anisotropy energy is phenomenologically thought to de¬ 
pend on the value of the c/a ratio, which is more or 
less equal for most of the compounds investigated in the 
present study. Therefore, the underlying mechanism is 
understood as a band-filling effect, which affects the spin 
orbit coupling (SOC) symmetry. This interpretation can 
be directly taken from Figs. 12 and 13. Increasing the 
SOC strength by varying the Y element through the 3d, 
4d, and 5d series increases the MCA energy by a factor of 
approximately 3 for Y = Fe, Ru, Os. The effect of band 
filling is deciphered by sweeping the Y elements along a 
series. Going from left to right in any set of compounds 
Mn 2 V^ xd -*Ga, the MCA is altered from preferred out-of¬ 
plane to in-plane orientation. The same situation holds 
for M^V*- 5 ^ Z compounds, whose preferred orientation 
is graphed in Fig. 13. Two compounds in this figure de¬ 
serve special attention, which are M^Ptln and M^IrSn. 
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Previously they were predicted to possess a non-collinear 
magnetic order just as M^RhSn, which was described 
in great detail in Ref. [46]. Although we did not con¬ 
sider non-collinear order in the present case, it might be 
worth while to point out that it is the Mn-atom on site 4 d 
that reacts by canting to the incipient spin-reorientation 
that is seen clearly in Figs. 12. Similar physics is ob¬ 
served in the famous Rare-Earth magnets Nd 2 Fei 4 B and 
Er 2 Fei 4 B. Manganese thus shares properties with Rare 
Earths [55-57], 


Group Number (y) 

3 4 5 6 7 8 9 10 11 



FIG. 12. (Color online) Calculated MCA energy of Mn 2 YGa. 
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> 

CD 

E 1 

a 

LU Q 


L2, 

1 

1 


Xa 
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1 
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1 
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1 


Mn2 ' ,ln -1- Mn 2 K ,5d, Sn 
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1 
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t 


1 

B 
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1 
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Mn 2 lrSn 


1 


H| Mn 2 PtSn 




Mn 2 Ptln H 


34.81 
23.21 
11.60 
0 5 

i— 

-11.60 


22 23 24 25 26 27 28 

Number of valence electrons N v 


29 


FIG. 13. (Color online) Calculated MCA energy of some cho¬ 
sen M 112 Y Z alloys: MnsGe, Mn 2 lrln, Mn 2 lrSn, Mii 2 PtIn, and 
Mii 2 PtSn. 


magnetic order does not change. A decrease or change 
of sign of the coupling constant JMn( 4 d)-Mn( 4 d) is typical 
owing to preferential anti-parallel coupling as caused for 
short Mn-Mn distances. Instead, the cubic to tetragonal 
transition even results in increased values, thus leading 
to an increase of Tq. Prima facie, the overall trend of 
increased Tq s cannot be traced back to a common mech¬ 
anism. In contrast, there is no exchange interaction that 
exhibits a similar behavior over the series. For instance, 
the findings in the case of Mn 2 FeGa are related to the re¬ 
duction of magnetic frustration, which is due to the com¬ 
peting anti-parallel interaction Mn(45)-Mn(4d), and par¬ 
allel interactions of Fe(4c) with both Mn(4d) and Mn(46) 
neighbors. Upon the tetragonal distortion a spin-flip of 
the Fe(4c) local moment is observed, due to the change 
of the sign of Mn(46)-Fe(4c) interaction. Upon the dis¬ 
tortion, the strength of Mn(46)-Mn(4d) interaction is al¬ 
tered by approximately A </;Mn(4b)-Mn(4b) = 56 meV in 
Mn2FeGa. 

A similar, but smaller effect is found in M^CoGa, indi¬ 
cating the magnetic frustration that had been present in 
the cubic phase and the weakening of the exchange inter¬ 
action Y (4c)-Mn(46) (Y = Fe, Co) due to the tetragonal 
distortion. This might be one of the contributions pro¬ 
hibiting the shape memory effect in M^FeGa. The main 
contribution to Tq, the Mn(4d)-Mn(46) exchange, does 
not suffer from the structural transition. Similarly, the 
influence of the distortion on the exchange causes the 
smaller contributions to be altered and some frustration 
to be diminished. For example, the Mn(4d)-Mn(4d) in¬ 
teraction vanishes, with preferred anti-parallel alignment 
in the cubic case, whereas the major Mn(4d)-Mn(46) in¬ 
teraction remains unchanged. 

Exceptions to the general observation of increased Tq in 
tetragonally distorted phases are Mn 2 NiGa, Mn 2 PdGa, 
and Mn 2 PtGa systems, in which the Xcs are reduced 
upon the tetragonal distortion. The significant reduc¬ 
tion (Fig. 14) is caused by a weakened Mn(4d)-Mn(46) 
interaction (Fig. 15) that may indicate an unstable mag¬ 
netic groundstate. A relation to the Heusler compound 
Mi^RhSn [46] can theoretically be established as these 
materials possess the same number of valence electrons, 
Ny. M^RhSn has been shown to exhibit a non-collinear 
magnetic groundstate due to canting of the different man¬ 
ganese moments [46]. Non-collinear spin configurations 
have not been considered in this work, so that canting of 
the spins or spin spiral groundstates in these materials 
may be present. 


C. Exchange Coupling and Curie Temperatures 

The details of the calculations pertaining to Curie tem¬ 
peratures (Tc) have been described by us previously [15]. 
Even though the in-plane next-nearest-neighbor distance 
decreases due to the tetragonal distortion, the next- 
nearest neighbor Mn(8c)-Mn(8c) coupling is still positive 
for the L2i-derived tetragonal phases; thus, the overall 


V. SUMMARY 

Using total energy calculations within density func¬ 
tional theory, we investigated in detail the response to 
tetragonal distortions for a large set of cubic Heusler com¬ 
pounds, Mn 2 Y^ 3d,4d ’ 5d ' ) Ga, and some other chosen mate¬ 
rials. We were able to single out the systems that remain 
cubic from those that favor a tetragonal structure. The 
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Pc 

Pt 

Tc,c 

Tc,t 

ATc,t-c 

-Emca 

K u 

Mn2ScGa 

87 

35 


464 


0.616 

1.62 

Mn2TiGa 

83 


557 





Mn 2 VGa 

94 


587 





Mn2CrGa 

97 

50 

578 

970 

392 

0.779 

2.46 

MnaGa 

96 

56 

221 

610 

389 

0.906 

2.7 

Mn 2 FeGa 

95 

56 

601 

848 

247 

0.359 

1.16 

Mn2CoGa 

93 

60 

928 

1124 

196 

0.236 

0.77 

Mn2NiGa 

35 

42 

1005 

750 

-255 

0.193 

0.62 

Mn 2 CuGa 

53 


1491 





Mn2ZrGa 

82 


207 





Mn 2 NbGa 

98 


289 





Mn2MoGa 

85 

65 

140 

335 

196 

0.636 

1.91 

MmRuGa 

95 

1 

619 

1315 

696 

0.564 

1.68 

Mn2RhGa 

15 

59 

576 

1351 

776 

0.322 

0.95 

Mn2PdGa 

7 

46 

809 

335 

-473 

0.040 

0.11 

Mn2AgGa 

24 


1240 





MnaHfGa 

89 







MnaTaGa 

96 







Mn 2 WGa 

83 







MnaOsGa 

96 


273 

1075 

802 

3.270 

9.72 

MnalrGa 

74 

5 

411 

1122 

711 

2.388 

7.02 

MnaPtGa 

23 

51 

799 

326 

-472 

-0.293 

-0.84 

MnaAuGa 

8 

26 

1027 

897 

-130 

-0.731 

-1.94 


TABLE IV. The calculated Curie temperatures in Kelvin of 
tetragonal Tc.t and cubic Tc.c parent compounds (taken from 
Ref. 15). The changes due to the tetragonal transformation 
are listed as ATc,t-c. Emca represents the magnetocrys¬ 
talline anisotropy energies in meV per formula unit, whereas 
the anisotropy constant, K u , is given in 


Group Number (Y) 

4 5 6 7 8 9 10 11 



Number of valence electrons N v 


FIG. 14. (Color online) The calculated Curie temperatures 
of the Heusler compounds containing Ga. The results shown 
here are obtained in the Mean-Field Approximation (MFA) 
and highlight the consequence of the structural relaxation. 
Squares correspond to cubic compounds and rectangles to 
tetragonal compounds. 


Sc Ti V Cr Mn Fe Co Ni Cu 



Number of valence electrons N v 

FIG. 15. (Color online) The evaluated effective exchange in¬ 
teraction parameters are shown and compared for all cubic 
and tetragonal cases. Thus, the underlying mechanism of 
the increase of To due to the distortion is visualized. Circles 
correspond to cubic compounds and rectangles to tetragonal 
compounds. 



FIG. 16. (Color online) A two-dimensional projection of 
the nearest-neighbor coordination of Mn(46) is shown. The 
symmetry of the coordination changes while undergoing the 
martensitic transition, whereas the nearest-neighbor distance 
remains unchanged as the volume change is on the order of 
only 1-9%. 


details of the total energy as a function of the distortion 
were found to be similar for materials exhibiting the same 
number of valence electrons. The magnetizations of the 
tetragonal alloys were found to be shifted to smaller val¬ 
ues, which we could attribute to changes of the itinerant 
local moments. This led to characteristic modifications 
of the Slater-Pauling curve. By means of partial den¬ 
sities of states, the changes to the electronic structures 
revealed the microscopic origin of the observed trends. 
As compared to the cubic parent phases, a strengthening 
of the exchange interaction between neighboring sites was 
observed, which resulted in an increase of the Curie tem¬ 
perature. Focusing our attention on the magnetocrys¬ 
talline anisotropy we observed an interesting trend that 
describes a spin reorientation over our series of com- 
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pounds, furthermore very large anisotropies are found 
for tetragonal Heusler compounds containing heavy tran¬ 
sition metals accompanied by low magnetic moments, 
which indicates that these materials may be promising 
candidates for spin-transfer torque applications. 
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